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In this preceedings paper we report on a calculation of graphene's Landau levels in a 
magnetic field. Our calculations are based on a self-consistent Hartree-Fock approxima- 
tion for graphene's massless-Dirac continuum model. We find that because of graphene's 
chiral band structure interactions not only shift Landau-level energies, as in a non- 
relativistic electron gas, but also alter Landau level wavefunctions. We comment on the 
subtle continuum model regularization procedure necessary to correctly maintain the 
lattice-model's particle hole symmetry properties. 
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1. Introduction 

Graphene, a one-atom-thick two-dimensional crystal of carbon atoms arranged in 
a honeycomb lattice, is a gapless semiconductor with an unusual massl ess Dirac- 
fermion band structure that has long attracted theoretical interest The low- 
energy properties of graphene are characterized by quasiparticle dispersion ^ linear 
in momentum and by vanishing density-of-states at the neutral system Fermi energy. 
The band eigenstates can be considiered as sublattice-pseudospin spinors and have 
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a chiral property which qualitatively alters the way in which electron-electron inter- 
actions influence electronic properties. In particular electron-electron interactions 
lead to a logarithimic enhancement of the Fermi velocit y in d oped and undoped 
graphene related to a lack of screening at the Dirac point SEEI. 
In the presence of a magnetic field, graphene's electronic structure also changes in 
a nontrivial way when compared to the non-relativistic two dimensi ona l electron 
gas (2DEG) case, leading to the so-called half-quantized Hall effect i n w hich 
the plateau values of the Hall conductivity are given by a xy = 4(n + l/2)(e 2 //i). 
Plateau conductivity values are separated by 4e 2 //i because of the fourfold degen- 
eracy due to valley and spin. In this paper we analyze the effect of electron-electron 
interactions on graphene's LL spectrum. We show that because of the chiral nature 
of graphene's band structure, interactions not-only shift Landau level energies but 
also alter Landau-level wavefunctions. 



2. Self Consistent Hartree-Fock Approximation 

In this section we use the self-consistent Hartree-Fock approximation (SCHFA) to 
study the effect of electron-electron interactions on the LL spectium of graphene 
within the massless Dirac-fermion(MDF) model. The low-energy properties of 
graphene can be adequately described by a MDF model: 

Hp = va-p (1) 

where a l, s are Pauli matrices acting on graphene's psuedospin degrees of freedom, 
p is a two dimensional vector measured from the Dirac points. (As we will discuss 
later this model requires especially subtle ultraviolet regularizations in order to 
yield physically correct predictions.) In the presence of a unifrom magnetic field 
B = Bz applied in a direction perpendicular to the plane of the graphene sheet [1] is 
modified by p — > 7? = p — (e/c)A where A is the vector potential with B = V x A. 
Defining the usual raising and lowering LL operators and a, with = (Ib/ \/2K)-k, 
where l B = {hc/eB) 1 / 2 = 25.66/yOB[T]) nm is the magnetic length, we can identify 
a zero-energy eigenstate given by acj)Q = and finite-energy chiral eigenatates n 

and n = —n with eigenenergies e„ = sgn{n)\J 2hv p^ T. i n the Landau gauge 
A = (0, Bx,0) 7 the corresponding eigenvectors are 

where 4> n ,x is a Landau-gauge eigenstate of a non-relatvistic electron gas and X 
denotes the guiding center degree of freedom within a LL. Projecting the interacting 
Hamiltonian onto the LL basis gives: 



Me-e = - ^v q : p(-q)p(q) 



(3) 
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where v q — 2ire 2 /e\q\ is the two-dimensional Fourier transform of the Coulomb 
interaction and 

P(-<1)= E ci XTS c n , xlTls ,6{ qv + X-X'y^ x '-^.F^ n ,{q), (4) 

nXn'X' 

In Eq.(d|) c i nXTS /c n xTs are creation/annihilation operator for particles in LL n at 
guiding center X for valley r and spin s. For notational simplicity we have assumed 
l B = 1; we will however restore these length units in the final results. F^ n ,(q) 
is referred to as graphene's relativistic form factor and captures the orbital and 
sublattice pseudospin character of the LL orbitals ^(n / and n' 0) : 

F n,n'(q) = ^[F\n\,\n>\{q) + sgn(nn')F lnl „ 1 ^ n , l _ 1 (q)], (5) 

F\n\,\n'\{'f) ls the well known form factor for an ordinary 2DEGEE^in the presence 
of a perpendicular magnetic field: 

F da) - \ V , ^ ~ (6) 

and the form factor for the lowest LL is just F$ (q) = F(q) = er q / 4 . Finally 

Here i^, - ™ are the associated Laguerre polynomials. 

In the Hartree-Fock approximation the effective single-particle Hamiltonian de- 
pends on the density matrix. In the case of Landau- level systems, the density-matrix 
is usefully parameterized as follows: 

K'^sA^ = W E (4,x', T ' )S ^, x , T>s )5( % +X'-X).exp-^ (8) 

The Hartree-Fock theory Hamiltonian is expressed in terms of A as follows: @ 

1 



(n,X,r, S |W e _ e |n',X',r', S ') = E E 



^ E W.-MW (p)A^, s , (?) exp 



L 2 

p 



exp^ x ' + X - X') 



here we restore Z# and define 

Vn u n 2 ,n 3 ,nM = V q F* uni {q)F*^{-q) (10) 

For the purposes of this paper we assume that the translational invariance is not 
broken. In this case A is non-zero only for q = 0. It follows that 

4S(?1 = KrO(W?=oS(\ ni \ - |n 2 |). (11) 



a The details of this calculation arc similar to the one in Ref. 10. 
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The 8 function which sets mixing between states with different values of \n\ can be 
seen to follow from spatial isotropy in the continuum model. Because both positive 
and negative values of n occur in graphene this restriction does not forbid mixing 
of states with positive and negative n by interactions. The fact that Landau-level 
wavefunctions are altered by interactions in graphene is the main difference between 
relativistic and non-relativistic cases. 

Assuming no broken translational symmetry, the first term in ([9]) is just the 
constant Hartree (electrostatic) potential which can be absorbed in the zero of 
energy. Dropping this term yields 

(n,X,r,s\H\n',X',r', S ')=- £ ^E<,»'®C(-«") A n£(0)^^') 

l»i|=l«2| q 

(12) 

The mixing between n and n LL is due to the fact that the self-energy is not diagonal 
in the pseudospin as can be seen by examining spatial isotropy consequences more 
closely. Isotropy is encoded in the form factors in the above equation: 

F* n ,(q)F* m (-q) = F* n ,(q)F* n (q) ~ cxp-Wl-M) (13) 

yielding an angular integral that is proportinal to S(\n'\ — \n\) not 5(n' — n). The 
other delta function gives X — X' . Assuming no spin or valley broken symmetries 
we can also assume that different valleys can be considered independently. This 
yields (suppressing the spin and valley indices) : 

(n\H e - e \n') = - J2 ^^v q F^A<l)F^ ni {-q)^^. (14) 
l"i|=|™2| q 

Notice that this matrix element is non-zero only if \n\ = \n'\. 

To compute the matrix elements we have to employ a LL index cutoff reminisent 
of the high-energy cutoff used in for MDF description of graphene at zero magnetic 
field. This ultraviolet cutoff plays a role because of the unbounded negative energy 
sea of the massless Dirac model. Following the procedure used at zero magnetic 
field, we choose a LL cutoff, a maximum value for |n|, M based on the physically 
natural cutoff of momentum at an inverse lattice constant scale, k c ~ 1/a, and on 
the semi-classical relationship between momentum and Landau-level index. This 
yields 

hv F k c = V2^VM. (15) 

Using k c ~ 1/a where a = 0.246nm is graphene honeycomb lattice constant we get 
a magnetic field dependent cutoff 

5000 

M (16 

B[T] v ; 

where B is the magnitude of the magnetic field. We can write the mean field Hartree- 
Fock hamiltoninan in the (n, n)-sector as a two level system: 

H MF =£+B-a (17) 
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Fig. 1. Hartree-Fock eigenvalues for a neutral graphene sheet in the units \/2hv/lB with a gr = 
2.2. In this calculation the energy of the n = Landau level is chosen as the zero of energy so that 
all n ^ eigenstatesoccur in particle-hole symmetric positive and negative energy pairs. 



where the pseudospin electric £ and magnetic B field in the (n, n)-space £ — 
-a gr {^/2){V n + V n )/2, B x = -<v(V?F/2)V„ and B z = y/n - a gr (^/2)((V n - 
V n )/2) , depend on the interaction matrix elements 



/b \2 \2 Ib \2 

Vn = -o \ ~(n\He-e\n), V n = -r\ -(n\H e - e \n), V n = -(n\H e - e \n), 
e z V 7r e V 7r e V 7r 

(18) 

a gr = e 2 /(Hvf) is graphcnc's coupling constant. 



3. Results and Discussions 

We must now address the regularization procedures that are necessary to extract 
useful predictive results from these Hartree-Fock calculations. The exchange energies 
V n and V n in Eq. (fT8|) both diverge with cutoff M like \f~M while B z diverges like 
InM and B x goes to zero like l/\fM. These large exchange energies are indeed 
physical because graphene's 7r-electron system has a high density of electrons, close 
to one per honeycomb lattice unit cell. This large energy is however neither easily 
measureable or of any great interest. Instead we want to use our Hartree-Fock theory 
to calculate the spacing of energy levels near the Dirac point and their dependences 
on the total Landau level filling factor. Progress can be made by simply choosing 
a convenient zero of energy as we do in the zero-field case. We propose using the 
energy of the n = Landau level, (0|i? e -e|0) evaluated for a neutral graphene sheet 
as the zero of energy. In zero field the analogous choice solves all problems, but that 
is not true in the massless Dirac model case as we now discuss. 
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In order to explain the problem which remains and our resolution of this prob- 
lem we go to a more microscopic level by considering properties not of the mass- 
less Dirac model but of the one-band nearest neighbour tight-binding model for 
graphene's 7r-orbitals. This model has the particle-hole symmetry property that 
when the Hamiltonian acts on a wavefunction which is restricted to one sublat- 
tice of the bipartite honeycomb lattice, it produces a wavefunction confined to the 
other sublattice. From this property it readily follows that eigenstates of the band 
Hamiltonian occur in positive and negative energy pairs which have opposite in- 
tersublattice phases and, importantly for the Hatree-Fock calculations, that the 
density-matrix of a neutral graphene sheet is just half of the trivial density matrix 
of a state in which all 7r-orbital states are full, 

K>7 tral = \ 8i>,i>. (19) 

This property is preserved in a magnetic field and implies that the role of generalized 
Hubbard model interactions at the neutrality point in Hartree-Fock theory is simply 
to shift all energy levels by an irrelevant constant. This property is independent 
of the dependence of the interaction on site-separation. When translated to the 
continuum model, this property implies that for the case of a neutral graphene sheet 
both £ in Eq.([TT|) should be independent of n and equal to (0\H e - e \0) and that B x 
should vanish. Although the regularization procedure discussed above recovers this 
result with errors that vanish with cutoff like \j\f~M, the particle-hole symmetry 
property is so essential to the observed properties of graphene sheets that these 
errors are uncomfortably large at practical values of M. We therefore propose the 
following regularization procedure for Hartree-Fock Landau level calculations for the 
massless Dirac model of graphene: i) Solve the Hartree-Fock equations for neutral 
graphene by setting Ao,o = 1/2, and A„ iM = (1 — sgn(n))/2 for n ^ 0. ii) For 
the neutral case choose (0|ff e _ e |0) as the zero of energy and set £ and B x to zero 
to compensate for the violation of particle-hole symmetry caused by a finite M 
cut-off iii) In the case of charged graphene sheets A„/ jn must be determined by 
a self consistent calculation in which £ and B x are replaced by the the difference 
between their neutral and charged system values. The Hartree-Fock energy levels 
of a neutral graphene sheet obtained by following i) are illustrated in Fig.([3|). The 
logarithmic dependence of all levels on cut-off M in this figure is expected to appear 
experimentally as a weak logarithmic correction to the \[B dependence of Landau 
level energies on field expected for non- interacting electrons. 
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